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Gravitomagnetic charge that can also be referred to as the dual mass or magnetic mass is the 
topological charge in gravity theory. A gravitomagnetic monopole at rest can produce a stationary 
gravitomagnetic field. Due to the topological nature of gravitomagnetic charge, the metric of space- 
time where the gravitomagnetic matter is present will be nonanalytic. In this paper both the dual 
curvature tensors (which can characterize the dynamics of gravitational charge/monopoles) and the 
antisymmetric gravitational field equation of gravitomagnetic matter are presented. We consider 
and discuss the mathematical formulation and physical properties of the dual curvature tensors 
. . . and scalar, antisymmetric source tensors, dual spin connection (including the low-motion weak-field 

' approximation), dual vierbein field as well as dual current densities of gravitomagnetic charge. It 

I is shown that the dynamics of gravitomagnetic charge can be founded within the framework of the 

■ above dual quantities. In addition, the dual relationship between the dynamical theories of grav- 

' itomagnetic charge (dual mass) and gravitoelectric charge (mass) is also taken into account in the 

M , present paper. 

, 
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] The gravitomagnetic monopole possesses the dual charge to the mass. In this sense, the gravitomagnetic charge 

^ ' can also be referred to as the dual mass [1,2]. Some authors referred to it as "magnetic mass" or magnetic-type mass 

I {magnetic-like mass) [3]. Historically, Newman et al. discovered a stationary and spherically symmetric solution (now 

^ known as the NUT solution) that contains a second parameter l besides the mass to, and regarded it as the cmpty- 

, space generalization of the Schwarzschild solution [4] . Demiansky and Newman found that the NUT space is in fact the 



I. INTRODUCTION 



^ • spacetime produced by a mass which has a gravitomagnetic charge [5] . Dowker and Roche independently rediscovered 
k>( \ this interpretation of the NUT parameter l [6]. The cylindrically symmetric solution of gravitomagnetic charge {i.e., 

^ ■ the cylindrical analogue of NUT space, which is the spacetime of a line gravitomagnetic monopole) [7] was considered 
' by Nouri-Zonoz in 1997. The various properties of NUT space attracts attention of some authors [3,8,9]. Lynden- 
Bell and Nouri-Zonoz gave a new and more elementary derivation of NUT space based on the spherical symmetry 
of gravitoelectric and gravitomagnetic fields of NUT space and the spherical symmetry of its spatial metric ^ab [8]. 
They proved that all geodesies of NUT space lie on spatial cones, and used this interesting theorem to determine 
the gravitational lensing properties of NUT space [8]. Miller studied the global properties of the Kerr-Taub-NUT 
metric [9] . Zimmerman and Shahir considered the geodesies in the NUT metric and found that the NUT geodesies are 
similar to the properties of trajectories for charged particles orbiting about a magnetic monopole [3]. Recently, the 
quantal field and differential geometric properties and related topics associated with the NUT space receive attention 
of several researchers. Bini et al. investigated a single master equation for the gauge- and tetrad- invariant first-order 
massless perturbation describing spin < 2 fields in the Kerr-Taub-NUT spacetime [10]. In their another paper, Bini et 
al. analyzed the parallel transport around closed circular orbits in the equatorial plane of the Taub-NUT spacetime 
to reveal the effect of the gravitomagnetic monopole parameter on circular holonomy transformations [11]. As far 
as the quantization of gravitomagnetic monopole is concerned, historically, many authors considered this problem by 
analogy with Dirac's argument for the quantization of magnetic charge [6,12,13]. In history, it was clear to see that 
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such an idea (as Dirac's theory) might not lead to a consistent quantum field theory for the gravitomagnetic charge 
[8] if the gravitomagnetic charge is truly present in the universe. This topic will be further taken into consideration 
in our other work regarding the quantum theory of gravitomagnetic monopole. 

Lynden-Bell and Nouri-Zonoz thought that it is natural to ask how general relativity must be modified to allow 
for gravitomagnetic monopole densities and currents [8]. They conjectured the generalization will be to spaces with 
unsymmetrical afRne connections which have nonzero torsion [8]. However, by analogy with the electrodynamics which 
is only changed to include d^j^T^^ = but does not come with an additional vector potential [8,14] if the magnetic 
charge truly exists, here we will suggest another alternative to this problem, i.e., constructing an antisymmetric dual 
Einstein tensor to describe the gravitational field produced by gravitomagnetic matter. Note that here we will not 
introduce extra unsymmetrical affine connections. 

This paper has a threefold purpose: to establish the theoretical basis (such as dual curvature tensors) on which 
the field equation of gravitomagnetic matter is suggested; to apply the gravitational field equation to some problems 
(such as the NUT metric); to consider other related topics (such as the dual metric, dual spin connection and dual 
current density). The physical reasons for the nature of the mathematical results obtained will be interpreted and 
discussed wherever possible. The paper presents a detailed and self-contained treatment of the fundamental problems 
of classical gravitomagnetic charge. 

The paper is organized as follows: in Sec. II, by considering the connection among the NUT solution, Einstein's 
equation and the dynamics of gravitomagnetic matter, it is shown that gravitomagnetic matter requires a gravitational 
field equation of its own; in Sec. Ill, we study the dual Riemann curvature tensors and dual Ricci curvature tensors. 
The relationship between these dual curvature tensors and the Einstein (and Ricci) tensor is also demonstrated. In 
Sec. IV, we take into consideration the comparison of the gravity theory (involving the gravitomagnetic charge) 
with the electrodynamics, which will provide clue to us on how to construct the gravitational field equation of dual 
mass (and hence its dynamics). Several fundamental subjects regarding the dynamics of gravitomagnetic matter are 
discussed in Sec. V, which are as follows: (i) by using the variational principle, the antisymmetric dual Einstein tensor 
is obtained with the dual curvature scalar; (ii) the antisymmetric source tensor of dual matter is constructed in terms 
of the four-dimensional velocity and the corresponding covariant derivatives; (iii) the dual spin connection which will 
arise in the field equation of dual mass is briefly suggested. The dual vierbein field, dual metric and dual current 
density of gravitomagnetic matter are taken into account in Sec. VI. The importance of these quantities derives from 
the fact that many dynamical variables (such as the source tensors) of physical interest regarding the gravitomagnetic 
matter can be expressed in terms of them. To consider the connection between the gravitomagnetic charge and the 
magnetic charge, in Sec. VII, we study the five-dimensional field equation of gravitomagnetic charge. In Sec. VIII, 
we discuss briefly the duality relationship between the dynamical theories of gravitomagnetic charge (dual mass) and 
gravitoelectric charge (mass). In Sec. IX, we conclude with some remarks. 

II. WANTING A GRAVITATIONAL FIELD EQUATION NECESSARY 

Is gravitomagnetic monopolcs required of its own gravitational field equation for dealing with the gravitational 
effects (properties, phenomena)? Before proceeding with the treatment of dynamics of gravitomagnetic matter, wc 
shall here consider a problem associated with the relationship among the field equation, the solutions and the Bianchi 
identities. In an attempt to investigate the metric resulting from the presence of the gravitomagnetic charge currents, 
we think that one may meet with problems concerning the choice of the "true" metric produced by gravitomagnetic 
charge from the nonanalytic solutions of Einstein's field equation. We argue that in history the gravitational field 
equation of gravitomagnetic charge might get less attention than it deserves. 

As is known to us all, in general relativity, the NUT metric (named after Newman, Tamborino and Unti) gni,{m, l) 
is the one describing the gravitational distribution produced by both the gravitoelectric charge (with mass m) and 
the gravitomagnetic monopole (with dual mass^ l) fixed at the origin of the coordinate system [5]. Since it is well 
known that the NUT solution 5^1/ (m = 0,l^0) satisfies Einstein's vacuum equation TZf_ii, = (and hence Einstein's 
equation G^i, = without matter), someone may hold that the dynamics of gravitomagnetic matter has already been 
embodied in Einstein'gravitational field equation, or that the gravitational field equation of gravitomagnetic charge 
has been involved in Einstein's equation, or that the gravitational distribution of both mass and dual mass is governed 
together by Einstein's equation. Thus from the point of view of these peoples it is concluded that there exists no 
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other gravitational field equation of gravitomagnctic matter than Einstein's equation, or that the potential suggestion 
of its gravitational field equation is not essential in investigating the dynamics of gravitomagnetic monopole (should 
such exist). However, even though it sounds reasonable, we think that these viewpoints may be not the true case. 
We believe that, in fact, the gravitomagnetic matter needs its own gravitational field equation rather than Einstein's 
equation, since the latter governs the spacetime of gravitoelectric matter (mass) only. In what follows it will just be 
concentrated on the field equation problem of gravitomagnetic charge by comparing it with the electromagnetic field 
equation of magnetic monopole. 

The NUT solution g^uim = 0,i ^ 0), which are the nonanalytic functions, truly agrees with Einstein's equation 
(empty-spacetime field equation), and it can truly describe the gravitational field near a gravitomagnetic monopole. 
This, however, does not means that all the nonanalytic solutions to Einstein's equation belong to the contribution 
of gravitomagnetic charge. Namely, Einstein's equation is overdetermined for determining the gravitational-potential 
solutions of gravitomagnetic charge. Moreover, it can be further concluded that the NUT solution gfj,v{m = 0, t 7^ 0) 
satisfying Einstein's equation is not necessarily the one only to Einstein's equation, just like the case that a solution 
satisfying a Bianchi identity is not the one only to this Bianchi identity. 

In order to clarify the above viewpoint further, let us consider an illustrative example. Since there are some 
analogies between general relativity and electromagnetism, and most of the essential features of dynamics of electric 
and magnetic charges share similarity to the gravity theories of gravitoelectric and gravitomagnetic charges, one can 
prove the previous interpretation to be sound by analogy with the Maxwellian equation in electrodynamics. It is 
well known that the electromagnetic field equation of electric charge and magnetic charge are written in the form 
d^T^"" = and 5^ J"^'' = j^^, respectively. Here, the dual electromagnetic field tensor JP""^ = le'^""''^ J^af} with e'"'"''^ 
being the Levi-Civita tensor^. It is readily verified that the stationary monopole solution (say, the Wu-Yang solution 
[15]"^) also agrees with the field equation of electric charge (Q^JF'"^ = withoTit the electric current densities). But 
this does not means that the Wu-Yang solution should necessarily the one only to = 0, since the equation 

{dij,T'^" = 0) is overdetermined for determining the monopole solution (for instance, it cannot set the magnetic charge g 
in the Wu-Yang solution). Likewise, it is apparently seen that the electric charge solution of d^T^^ = automatically 
agrees with the electromagnetic field equation of magnetic charge (i.e., d^^^^ = without magnetic charge densities 
and currents). But it is clearly not suitable to say that the equation d/^f^^" = is also the electromagnetic field 
equation of electric charge current . In fact, not all the analytical solutions which agree with df^T'"'' = satisfy 
the equation d^J^'^'^ = j^. The reason for this may be as follows: for the latter field equation, the former equation is 
merely a Bianchi identity. Similarly, for the magnetic charge, d^J^i^" = should also be viewed as a Bianchi identity, 
which may also be seen as follows: if the dual electromagnetic field tensor is defined to be !F^'^ = d^A^ — d^A^, 
and the electromagnetic field tensor Tf^" expressed in terms of the dual electromagnetic field tensor takes the form 
J^f^^ = —^€'^'^°'^!Faf3, then the electromagnetic field equation of electric charge {i.e., d^!F'^" = without the electric 
current densities) can just be rewritten as the form of a Bianchi identity. Thus, we show that the equation d^^J^^'^ = 
in the framework of dual 4-D potential vector A'^ is truly a Bianchi identity. Since d^J^^'^ = is only a Bianchi identity 
for the magnetic charge, it cannot be the electromagnetic field equation of magnetic monopole (similarly, Einstein's 
equation in empty spacetime is not the gravitational field equation of gravitomagnetic matter, cither). This, therefore, 
means that not all of the solutions satisfying this Bianchi identity may also satisfy d^^'^^ = j^, i.e., not all these 
solutions may surely belong to the effect of magnetic monopole (similarly, not all the empty-space generalizations of 
the solutions of Einstein's equation belong to gravitomagnetic charge, i.e., Einstein's field equation cannot govern 
independently the gravitational distribution produced by gravitomagnetic matter). In a word, gravitomagnetic matter 
should have its own gravitational field equation. 

In view of the importance of gravitational field equation of gravitomagnetic monopoles, it is certainly desirable to 
construct a theoretical framework, within which the field equation (and hence its dynamics) can be established. First 
and foremost, we should study the dual curvature tensors in Riemann space. 



^Note that, in this paper, for tlic electrodynamics, the Levi-Civita tensor refers to that in Minkowski flat spacetime, while for 
the gravity theory, e^"™'' stands for the Levi-Civita tensor in the curved spacetime. 

^The Wu-Yang solution is of the form (in spherical coordinate system) Ar = 0, Ag = 0, A^f, = which is a nonanalytic 

function. 
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III. DUAL CURVATURE TENSORS IN RIEMANN SPACE 



In Sec. II, we showed that the investigation of gravitational distribution in the vicinity of a gravitomagnetic 
chargc/monopolc should start with the needs of a gravitational field equation. In order to discover the gravitational 
field equation of gravitomagnetic matter, we will proceed in a manner entirely analogous to what has been done for 
that of magnetic monopole (i.e., first we define the dual electromagnetic field tensor J^i^" , and then we obtain the 
quantity df^T^'^ from the variational principle with the dual Lagrangian density). 

To lay a mathematical foundation for the field equation of dual mass that we will develop in the next section, first 
we consider the properties of dual curvature tensors in Riemann space. The right- and left- dual Riemann curvature 
tensors may be defined as follows 

n*^r^, = ^e^/^n^rxa (right dual), *7^^,^, = ^e^r^^TZxau^u (left dual), (3.1) 

and consequently, the corresponding dual Ricci tensors take the form 

= 9""Kru>. = 'l^.'^^T^f^rXa (right dual), =5™*^^™. = ^e^^""7eAac.. (left dual), (3.2) 

respectively. Note that ^^'^^'^'R-^rXa = 0, i.e., the right-dual Ricci tensor TZ*^^ is identically zero. So, we think that 

7?.*,^ may have no important value in dealing with the dynamics of dual mass and therefore we abandon if* and only 
adopt the left-dual curvature tensors in the following discussion. In particular, in what follows, the left-dual Riemann 
tensor and its corresponding dual Ricci tensor will be denoted respectively by TZfiTwv and TZ^^v, i.e., 



1 ^ 1 

2^ 



It should be pointed out that in this section, without loss of generality, we shall restrict ourselves to the left- and 
right-dual curvature tensors only. For other dual curvature tensors (such as mid-dual and two-side dual tensors) one 
may be referred to the appendix to this paper. 

To consider further the properties of dual curvature tensor, we will present the following three interesting and useful 
identities regarding the above-mentioned dual curvature tensors: 



Here TZ'"'^ = TZ"^^ ^, TZ'^^ = Ti^^ Xu- The proof of the first identity in Eq.(3.4) is given in Appendix 2 to this paper. 
Now let us look at the second identity in Eq.(3.4), which can be easily verified as follows: 

1 ^ „ 1 



where for the calculation of e^'^^^t^J^'^ , readers may be referred to Appendix 1. In Appendix 1 the Ricci tensor is so 
defined that 7^''^^^ = TZ^^. By using TZ^"" ^ ^ = VJ^^ , the second identity in Eq.(3.4) is proven correct. With the help 



of the formulae^ in Appendix 1, one can arrive at 



*It should be pointed out that in the published paper [Gen. Relativ. Gravit. 34, 1423-1435 (2002)], however, we adopted the 
right-dual form rather than the left-dual one. Such a situation was caused by the fact that the form of the Riemann curvature 
tensor TZ^i,a0 applied in the GRG paper (see, for example, Eq.(7) in the above-mentioned GRG paper) refers to the form TZag^iv 
expressed in the present paper. 

^It follows from the formula presented in Appendix 1 to this paper that the inner product of two Levi-Civita ten- 
sors/symbols (where the summation is to be carried out over the repeated index /u = 0,1,2,3) reads g'^''"'^^^ Actij _ 

- {g^'^g^'-g^^ + g^^g^'^g'^'^ + g^^g^^g"" - g'-^g^^g^" - g^'g^'^g^'^ - g^^g^^g^'^) ■ 
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Thus we obtain 



2 ~ 4 



It follows that 



(3.7) 



(3.8) 



which is the third identity in Eq.(3.4). Furthermore, by the aid of the formula in Appendix 1, one can show that the 
fourth identity in Eq.(3.4) holds also. 

In the following we will briefly discuss the difference between 71^,^ and TV For the Riemann tensor, even though 
TZaiivf} = —T^oitifiv holds, the relation TZafivis = —T^nav/) is no longer valid due to the nonanalytic property of the 
metric tensors, which can be seen from the following calculation: 



(3.9) 



Thus we have TZa^^/s = TZ^affu + {dud/3 - d^d^) 5a/n- With the help of 7?.^^ = g"^'Ra,j.ui3 and 7^'^^ = 5"^^^pa/3i/, one 
can arrive at 



(3.10) 



which is the relation between TZf^i, and TZ' ^i,. Apparently, if is an analytic function, then we have TZf^, = TZ' f^,. 
It should be noted that in the dynamics of gravitomagnetic matter, the metric is often nonanalytic. For this reason, 
one should distinguish TZ'^i, from TZ^i, in the treatment of the problems where both mass and dual mass are involved. 



IV. THE CONNECTION BETWEEN THE GRAVITY THEORY AND THE ELECTRODYNAMICS 

It is useful to consider the analogies between the gravity theory and the electrodynamics in developing the grav- 
itational field equation of dual matter. This deep connection between them will make the establishment of gravity 
theory of gravitomagnetic matter more enlightening and possible. 

The first analogy between them is that the Lagrangian density of electromagnetic fields can be constructed in terms 
of the dual field tensor J^/^v, i-e., 

-Jjr^.F"- = - (^-ij-^..^^'^^ with = ie^,"'='.F„^, (4.1) 

the gravitational analogy to which arises also. By making use of the identities of Eq.(3.4), it is shown that 

^e^'^'^-'TZr^,^ = -TZ (4.2) 

with TZ being the curvature scalar. 

The second analogy is that one can obtain the term d^^^" that appears on the left-handed side of the electro- 
magnetic field equation of magnetic charge currents from the variational principle with the dual Lagrangian density 
C = —jT^v^^'^- Note that here the variational principle is applied to C with respect to . In the similar manner, 
for the gravity theory, the dual Einstein tensor G^v can also be derived via the variational principle with the dual 
curvature scalar TZ, i.e., 

S I V^TZdn = ^ G^^ = TZi,^ - n^u,, (4.3) 

JQ. 
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where ft denotes the spacetimc region of volume integral and dfl is the volume element. The dual curvature scalar 
n is defined to be = g^''"R-f,^ = ^e'^^^^'TZxarix with TZ^^ = ^e^^'^'^TZxcrT,.. The relation (4.3) will be developed in 
more details in Sec. V. 

The third analogy between the gravity theory and the electrodynamics may be seen from the following three set of 
expressions. The relation between the electromagnetic field tensor jF^i^ and its dual JF^^ is 

In the meanwhile, there exist the similar relations between curvature tensors and their dual. For example, 

and 

= ^e/"^7^A.r., \ {n^.+Ti'^^) = " ^""^A<xr.- (4.6) 

Note that the first expression of Eqs.(4.6) is due to the definition of the dual Ricci tensor TZ^^i, in (3.3), and the second 
expression of Eqs.(4.6) can be derived from the second identity in Eqs.(3.4). 

We have, therefore, investigated in detail the properties of dual curvature tensors in Riemann space. Thus, the 
foundation of the dynamics of gravitomagnetic matter is laid in the preceding sections. Moreover, we have also 
discussed the problem of the requirement for a gravitational field equation of gravitomagnetic matter of its own 
in studying both the gravitational field distribution produced by the gravitomagnetic charge and the equation of 
motion of gravitomagnetic monopoles. In the following section, we will derive the gravitational field equation of 
gravitomagnetic matter from the variational principle with the dual curvature scalar. 



V. GRAVITATIONAL FIELD EQUATION OF GRAVITOMAGNETIC MATTER 

What in gravity theory may be caused by the existence of the gravitomagnetic monopoles? Will it lead to unsym- 
mctrical metric and/or to modify Einstein's field equation? As the topological charge in spacetime, gravitomagnetic 
charge will inevitably result in the nonanalytic parts in the metric functions. In this sense, gravitomagnetic monopole 
does not introduce the so-called extra gravitational potentials. So, here we need not employ the unsymmetrical met- 
ric in the treatment of the dynamics of gravitomagnetic matter. In order to determine the nonanalytic parts of the 
metric functions caused by the gravitomagnetic monopoles, we should first suggest a gravitational field equation of its 
own. A systematic procedure for the gravitational field equation of gravitomagnetic matter will be developed in this 
section. It will lay a physical foundation for the mathematical formalism in Sec. IV. If the gravitomagnetic charge 
is not present, then the contribution of this gravitational field equation without source term in general relativity is 
easily shown to vanish, on account of its status of Bianchi identity. But, it will play an essential role in treating the 
gravitational distribution problem of gravitomagnetic matter (should such exist). For examples, we shall illustrate 
the application of this field equation to the NUT solution and the equation of motion of gravitomagnetic monopoles. 



A. Dual curvature scalar and variational principle 



With the introduction and motivation in the previous sections, we are now in a position to find a gravitational field 
equation of gravitomagnetic matter. For this aim, one should first construct a Lagrangian density that can describe 
the gravitational field with nonanalytic property. Such a Lagrangian density should be a one that will be vanishing if 
the metric function is analytic. The only one we can choose is just the dual curvature scalar Ti. So, the dual action 
I' of the gravitational field produced by the gravitomagnetic matter may be of the form [1,2] 



Jn Jn L 



(5.1) 



which is a four-dimensional volume integral of the dual curvature scalar TZ. Note that the dual Ricci tensor 72.^^ is 
asymmetric in fi, v. So in Eq.(5.1), we introduce a certain parameter ( to rewrite the dual action /'. The variation of 
the action I' is 
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SI' = [ 6 iV^g'^'') [(1 - 0%. + CR-u^ 
The first term on the right-handed side of Eq.(5.2) can be rewritten as 



dfi + 


/ ^/^gg^'' 




Jn 



(1 - C)5^^. + c^^ 



dO, 



(5.2) 



siV^g'"') (1 - C)^^. + C^^^M 



dO 



-g 



(1 - C) ( 7^^. - IgM + C (t^.a. - Ig.f.'^ 



= J V-g - ^giJ.iyT^ - C (T^iJ.<y - 



(5.3) 



Now in the following we will calculate the second term on the right-handed side of Eq.(5.2). For convenience, we set 
A = \/— 55'"' (1 ~ O^T^nv + Q^T^vix dn. The infinitesimal variation of the dual Ricci tensor STZ^v is given 



(5.4) 



In Appendix 3, we show that the variation \/~gg^^\^fi ^"^^ STZxcrTu<i^ if' actually a surface integral over the boundary 
of the volume, namely, it can be transformed into a four- dimensional volume integral of a divergence. Since the varia- 
tional principle assumes that the variation 6g^" vanishes on the boundary, it follows that J^^ ^/—gg'^^^^fj, ^'^'^ STZxari^dfl 
has no effect on the variation of the action and therefore gives no contribution to the derivation of the dual Einstein 
tensor G^,/. We can, therefore, ignore this term^. Thus, the term STZ^vd^ in A can be reduced to 



JO ^ Jn 

By using the formula in Appendix 4, Eq.(5.5) can be rewritten as 

= {-gle^^'^^Uxaru + Igaffe'^^'^^nxaru^ Sg^^d^l 

In the same fashion, the term (/(/^''JT^.^^dO in A can be rewritten in the form 

/ ^g'^'^sn^^dn = - [ ^ (n^, - Ig^^a) Sg'^'dn = - / 

Jn Jn \ J Jn 

By introducing two parameters ^ and ^, Eq.(5.6) and (5.7) can be rewritten as 

^ V^g'^-'siz^^dfi = -(1 - ^ (n^. - Igt^.n^ Sg^"" dn ~^J^V^ (^.^^ - ^s-a*^ 

and 



(5.5) 



g { - ^gu^n ) ^^^'do. 



(5.6) 



(5.7) 



(5.8) 



^ ^/^5''^(5^,^d^^ = -{!-<;) j^^g (u,^ - ^g^^ll^ dg'^'^dO. (^/^^ - \g^^^'^ dg^^'^dn. (5.9) 



^Someone may argue that because of the symmetric property of the ChristofFel symbol in the indices /i, v, it is apparent to 

see that the integrand ^/— Sfl'"'e^ ^"^ STZxctti^ is strictly identical to zero, and it may drop out of the variation 51' immediately. 
For the dual formalism, however, this viewpoint is not appropriate for determining the dual Einstein tensor. To derive 
the dual Einstein tensor in a stringent way, if \/— SS^'^e/j ^"^ STZxarvAQ, can be truly ignored, we should first show that 
^J—gg'^'^t^ ^"'^ &R,\aTv can be rewritten as a covariant divergence of a certain vector indeed. 
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Hence, by using Eq.(5.8) and (5.9), one can express the second term on the right-handed side of Eq.(5.2) (i.e., 

^ = In V^g""" [(1 - Qsn^,^ + csn^^] dn) as 



(1 - c)6n^. + C6Ti.^ dn = [-{!- 0(1 - - a] / t^^.. - ^g^.T^ ] Sg^"dn 



+ [-(1 - C)^ - C(l - ?)] / 7^.^ - 25-/^^ ) ^ff^'dfi. 



Namely, we have 

A = -j^^ (n^^ - ^g^^n^ 6g^"'dn + [C + ^ - C(^ + ^)] {n^. - ^.^) Sg'^'^d^. 

It follows from Eq.(5.2), (5.3) and (5.11) that the variation of /' is finally given 

51' = K - C(C + ';)]j^V^ - ^.m) Sg^^dn. 
If we set / = -^zr^p^I' , the expression for the variation of / reads 



51 



-gGi^Jg^^dn with G^^ = 11^,^ -U^^,, 



(5.10) 



(5.11) 



(5.12) 



(5.13) 



where G^i/ may be referred to as the dual Einstein tensor, which is the resulting one for the given dual action /. 



In what follows we will simplify the form of the dual Ricci tensor 7^^^ ~ ^^^^'^^T^Xaru, which yields 



1 d ( dg^x dgav 



4 ^ dx'' \ dx'^ dx^ 



(5.14) 



Here any symmetric parts in TZxarv have simply dropped out of 7?.^,^ because is completely antisymmetric in 

A, (7, T. It is physically interesting that the (/iO) component of the dual Ricci tensor is TZ 
which is somewhat analogous to its electromagnetic counterpart jf^F'^ 



— If Actt d ( dgox dg^g 
^'O 4"^ dx^ \ dx" dx^ 

£^ {dxA„ - d^Ax). This, there- 



fore, means that G^v may truly serve as a dual Einstein tensor and would appear on the left-handed side of the 

gravitational field equation of gravitomagnetic charge. Clearly, the formalism for the dual Einstein tensor which we 
have developed here is rigorous. It may contribute to a better understanding of the dynamics of gravitomagnetic 
matter. 



B. Antisymmetrical source tensor 

The preceding subsection shows that the dual Einstein tensor G^i, is a second-rank antisymmetric one. So, the 
source tensor of gravitomagnetic matter, which appears on the right-handed side of the gravitational field equation of 
gravitomagnetic charge, should also be an antisymmetric one [1,2]. Such antisymmetric tensors which arc constructed 
in terms of the four-dimensional velocity vector Ui, and its covariant derivatives with respect to the spacetime coor- 
dinates may be only the following two forms 

IC^u = D^f/, - D,U^ = d^U, - d^U^, H^, = ^e^^"^lCc,0. (5.15) 

So, the second-rank antisymmetric source tensor of gravitomagnetic matter can take the following form {i.e., the 
linear combination of IC^i, and H^u) 

= KilC^„ + K2'Hn„, (5.16) 

where ki and K2 denote the two combination coefficients. 

For the Fermionic field tp, its four-vector velocity (or current density) is = e^i'tpjai', where denotes the 
vierbein field with /i, a = 0,1,2,3 and 7a 's stand for the Dirac matrices in the flat spacetime. The second-rank 
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antisymmetric source tensor of the gravitomagnetic Fermionic field ip can be constructed according to the above- 
mentioned definition^. The four-dimensional velocity (or momentum) of the Bosonic field is Ufi = {ip*dfj,(p — ifdntp* ), 
and the antisymmetric tensor /C^^ is therefore K.^^ = i {d^ip^d^^p — d^ip*dij,(p), which can be used to express the 
antisymmetric source tensor 5*^,^ of the gravitomagnetic Bosonic field. 

It should be noted that the source tensor of gravitomagnetic matter can also be expressed in terms of the dual 
current density (velocity) and the dual vierbein field e^", which will be discussed in Sec. VI. 

C. Gravitational field equation of gravitomagnetic matter 

In accordance with the above discussion, we discover the gravitational field equation of gravitomagnetic matter as 
follows (covariant form) 

Gfj.v = Sfj^^. (5-17) 

Note that Einstein' field equation is of symmetry in indices. In contrast, here the gravitational field equation is an 
antisymmetric one [1,2]. 

So far wc have obtained the gravitational field equation of gravitomagnetic charge of its own. Now let us discuss the 
equation of motion of a gravitomagnetic monopole as a test particle (or a particle system consisting of gravitomagnetic 
monopoles). According to the gravitational field equation (contravariant form) 

G^"" = (5.18) 

we set a second-rank antisymmetric tensor Z'^'^ = S^^ — G^^ , the covariant divergence of which is expressed by 

which means that Z^'^ = --^=C^^ with C^" being a certain constant tensor, or a tensor whose divergence vanishes. 
Because of Eq.(5.18), C^" = 0, i.e., Z^^" = 0. Note that even though the equation Z^^'^ = can be derived 
straightforwardly from Eq.(5.18), here wc would not think of it as the result of Eq.(5.18). On the contrary, wc prefer 
to consider it as the result of the covariant divergence of the field equation Eq.(5.18). This, therefore, implies that 
the equation Z^" = is inevitably in connection with the kinematic equation of a system of particles consisting of 
gravitomagnetic charges. The following calculation will confirm this interpretation: Si^^Uv in Z^^Uv is expressed by 

S'^^'U^ = KiU^ (D'^t/'^ - D^t/'') + l^ae^'^^'^t/. {daUp - df,Ua) 
nrr/i i 



ds 2 

= -Ki 



+ lK2e^''"'^U^ {daQpx - dpgax) + \K2e^'"'^U, {90xda - Qaxdp) U\ (5.20) 



It is noted that the second term ^K2e'^^°'^U^ {da9px — dpgax) on the right-handed side of (5.20) contains an 
expression ~ ^K2e'^''°'^ {dag/so — d/sgao) U^, which is just a gravitational (gravitomagnetic) Lorentz force, the elec- 
tromagnetic counterpart of which is ge^'^°'^ {daA/} — d^Aa) with g representing the magnetic charge. The fact 
that ~ \k2^^^'"^^ {dagpo — d/sgao) Uv resembles the expression for the Lorentz force in electrodynamics means that 
the equation Z'^" = contains the dynamical equation of motion of a particle system composed of gravitomagnetic 
charges. 

The antisymmetric field equation (5.18) of gravitomagnetic matter governs the gravitational field near the gravit- 
omagnetic monopoles. It does not introduce new extra gravitational potentials. The role of this equation presented 
here is to determine the nonanalytic parts of the metric which is caused by the existence of the gravitomagnetic 
monopoles. If the gravitomagnetic charge does not exist in the universe, the vacuum case of Eq.(5.18) is just the 
Bianchi identity of Einstein's field equation. 



'^The tensor /C^^^ of Fermionic field is given as i [(d^ip^ -y^ip — {di,^p) ^tii^] + # (Ti^^^m — liidi^) i> + ii> {d^-y„ — d^-y^) tp- In the 
papers [1,2], the third term, iip (9^7i/ — d^y^) ip, was unfortunately neglected. 
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D. About NUT metric 



The NUT metric is the one produced by the stationary gravitomagnetic monopole at the origin of the coordinate 
system. Here we will discuss the relation between the NUT parameter b and the gravitomagnetic charge. The 
contravariant-form field equation (5.18) of dual matter can be rewritten as 

_ ^ 'Ynv giiv _ 1 ^^nv _ fi^n^ ^5 21) 

For a single gravitomagnetic monopole, the strength of which is t', one may assume that T" q = l'6{x). The four- 
dimensional volume integral of the term on the right-handed side of the field equation TZ^ q = -^T^ q is 



/ 

Jn 



^ { ^Lt" ) dfi = l'. (5.22) 
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In what follows, we will consider the equation 72.° q = -^=L'd{x), where TZ^ q = jc'^^'^'^dr {dagox — dxgao) and e° 
_^g.0A<7r jjgj-e gOAar jg ^j^g Levi-Civita tensor in the flat Minkowski spacetime. It is readily verifled that this equation 
can be rewritten as 

V • Bg = 4t"5(x), (5.23) 

where t" = l'5{x'^), S{x) = d{x'^)S(x). Here the gravitomagnetic field strength is defined as Bg = V x g with the 
gravitomagnetic vector potentials being g = {goi, go2, goa)- 

It is well known that the gravitomagnetic vector potentials of NUT metric (without mass) can be expressed as 

g=(0,0,-2/V)^^^^:^) (5.24) 
\ r sm t' / 

in the spherical polar coordinate system, namely, go^ = —2f'^{r) ''^^~?°^f\ where the function /^(r) is /^(r) = 
1-2 (p^) • Thus it is easy to obtain 

Vxg= ^ep + -^^^^ ^-r^^o- 5.25 

r or rsmO 

It should be noted that the divergence V • (V x g) no longer vanishes at some spatial points (singularities) since here 
the gravitomagnetic vector potential is nonanalytic, which results from the presence of the gravitomagnetic monopoles. 
Hence the three-dimensional volume integral of V • (V x g) is 



X 



V • (V X g) dV = 1^ (V X g) • dS = -Sttl. (5.26) 
^ s 



It follows from Eq.(5.23) and (5.26) that the relation between c" and the NUT parameter (, is given 

l" = -2TTi. (5.27) 

It should be emphasized that the Schwarzschild solution also satisfies the field equation (5.18) without source term, 
i.e., = 0. Someone may thus argue that the Schwarzschild solution is just the empty-space generalization of the 
solutions of Eq.(5.18). This situation is similar to the case of the NUT solution g^v{fn = 0,t ^ 0), which was also 
thought of as the empty-space generalization of the solutions of Einstein's equation [4] . 

How can we look upon the two solutions 5^1/ (m = 0, i ^ 0) and g^vim ^ 0, t = 0) in the NUT spacetime? It is 
believed that the vacuum equation (G'*'' = 0) of dual matter may be viewed as the Bianchi identity of Einstein's field 
equation, and Einstein's vacuum equation is in turn the Bianchi identity of the field equation (5.18) of gravitomagnetic 
matter. For this reason, even though the solution 5^1/(771 7^ 0, t = 0) is truly the one to Einstein's field equation, it can 
also agree with the equation G'*'' = (Bianchi identity). So, it may be argued that the solution g^„{m ^ 0,l — 0) is 
the empty-space generalization of the solution of vacuum equation G'^'^ = 0. Likewise, the reason for Newman et al. 
to think of g^vim = 0, t 7^ 0) as the solution of Einstein's vacuum equation [4] lies in that Einstein's vacuum equation 
is just the Bianchi identity of Eq.(5.18) of gravitomagnetic matter. In fact, gixv{m = 0, i 7^ 0) is only the one of the 
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solutions of Eq.(5.18), and furthermore, all the solutions satisfying Eq.(5.18) will also agree with Einstein's vacuum 
equation (Bianchi identity). In this sense, the gravitomagnetic charge is truly required of its own gravitational field 
equation for us to treat its dynamical problems. 

In addition, we consider briefly the stationary cylindrically symmetric exact solution of field equation (5.23). Sup- 
pose that the form of linear element describing the stationary cylindrically symmetric gravitomagnetic field is given 

by 

ds^ = {dx°f - da;^ - dy^ - dz"^ + 2go^{y)dx°dx + 2gQy{x)dx°dy, (5.28) 

where we assume that the gravitomagnetic potentials go^ and goy are the functions with respect to y and x, respec- 
tively. For the cylindrically symmetric uniform gravitomagnetic field produced by the gravitomagnetic charges (with 
a nonvanishing surface density of um) present only in the x-y plane of 2 = 0, the gravitomagnetic field, Bg, which is 
defined as V x g, may be written 

(Bg), = ^. (5.29) 

It follows from Eq.(5.29) that the direction of Bg is parallel to the 2;-axis. The metric components of the uniform 
gravitomagnetic field, go^, 9oy, are therefore readily obtained as follows 



90x 



Ba 



-y, goy 



with Bg = ^. 



In order to obtain the contravariant metric g^^ , we calculate the inverse matrix of the metric {guv) 
is given as follows 



(5.30) 



and the result 



1 + flox + gty 



( 1 

.901/ 

V 



gvix 
-(l + 5o\) 

goxgoy 




goy 
goxgoy 





-(1 








■ gox 



\ 



+ gly)) 



(5.31) 



E. Dual spin connection 

Here we will consider the covariant derivative of Dirac field composed of gravitomagnetic charge in the local frame. 
Assume that an observer moves with proper three-acceleration a and proper three-rotation Co inside an inertial frame 
of reference. Note that here the quantities and symbols just refer to those in the reference [16]. It is believed that the 
covariant derivative of gravitomagnetic- charge Dirac field is given as follows 

D& = de.-^^a^'^e^.^'T^^,^. (5.32) 
Now we discuss the dual spin connection —i^<^^^ £ ^-'^^^ xf&- using the following relations [16] 



ijO ~ 1 _|_ aoc ' 010 ~ 100 ~ i i aoc ' fiPO ~ 000 ~ 

la'^ = ~, ia°^=a, (5.33) 



one can verify that 

"Tel rij0 + 2£7 £(,c ToiO 



16 V" """^ ' ' — * "'"y ]^ _|_ a x I Ig be ^ijfe ^ + ^bc ^2 

1 / 1 ^ ^ i 



uj ■ a ■ 



1 + ^ V4c 2nc2 



a • S , (5.34) 



where the Latin indices run over 1, 2, 3. Here a and S denote the velocity operator and spin operator of Dirac field, 
respectively. 

It is well known that a Dirac particle in a noninertial frame of reference will undergo the spin-rotation coupling 
(w-S) [17,18] and velocity-acceleration coupling (a- a) [16]. However, it follows from (5.34) that for a gravitomagnetic- 
charge Dirac field, it experiences velocity-rotation coupling {u ■ a) and spin-acceleration coupling (a • S), which are 
just the dual interactions of the above two couplings. 
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VI. DUAL METRIC AND DUAL CURRENT DENSITIES 



A. Dual vierbein field and dual metric 

It can be shown that in the weak gravitational field, Dirac matrices 7^'s is rewritten as 7*^ = (^1 + l^i 7* = 

+ 7m (i = 1, 2, 3) [1]. Here /3 and 7^ are the Dirac matrices in the flat Minkowski spacetime. (f) and can be 
viewed as the gravitational scalar potential and gravitomagnetic vector potentials. By using 7** = e'"*7a (a denotes 
the indices of coordinate in Minkowski spacetime) , one can obtain the relationship between the vierbein field e^" and 
the gravitational potential and gravitomagnetic vector potentials: e°° ~ 1 + ^ and ~ 3^° {jjL = 1,2,3). Here the 
dual vierbein field e^a is so defined that e^a and e^a satisfy the following relation 

d^iua - d„e^ta = ^e^L-"^ {dae^a - dfjeaa) ■ (6.1) 

Moreover, the dual Dirac matrices can be defined as 7^ = e^al""- Multiplying the two sides of the above equation by 
7", one can arrive at 

- = ^e^^."^ {dal/} - 9/37a) , (6-2) 

which is the relation between the dual Dirac matrices 'y^ and the regular Dirac matrices 7^. 

The connection between the metric 5^,^ and the flat metric r/ab is Qfiv = e^ej'riab- Correspondingly, the relation 
between the dual metric (jf^i, and the flat metric rjab may be = e^ej'rjab- It should be noted that even though we 
introduce the dual metric g^,^, in general, indices of tensors are raised and lowered still with the metric tensor g^^, 
rather than with the dual metric 5^1/ . 

In this subsection, we will discuss the dual covariant derivatives, V^V) of Fermionic fleld ip [19]. The dual spin 
connection of V^i^V = (J^ij. + ip is deflned by 



B^ = ^efsce\,^s^%y', (6.3) 



where a'"^ = | [7''. 7''] . It is essentially signiflcant to construct the Lagrangian density of Fermionic- fleld gravitomag- 
netic charge, which may be given 



I ^e^" [^7aV^^ - (V^V;) 7a V'] - m^iV-j , (6.4) 

where m denotes the dual mass of gravitomagnetic monopole, and S^^'ip = (j)^ — iB^^ 'ip. The corresponding classical 
field equations of tp and ■tp are as follows 



i 1 



=5^ iV^^e'^^'^a) - \ (V^V^) e^"7a - m^i = 0. (6.5) 



It follows from Eqs.(6.5) that the conservation law of the dual current density is 

1 



(V^i^g'^"7aV') = 0, (6.6) 

where the dual current density of dual mass is 

f = i^e^'^j^i^ = it^-j,, (6.7) 

which will be considered further in the following. 

In addition, the Lagrangian density of the Bosonic-fleld gravitomagnetic matter may be of the form 

jOb = [e""" {d^^*daV + daip*d^ip) - mVV] • (6-8) 
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The corresponding classical field equations are 



^ 5^(y^e'^''5„<p*)+mV*=0. (6.9) 



V-5 

It follows that the conservation law of the dual current density takes the form 

1 



[7=56"" {<P*da^ - <pda^*)] = 0. (6.10) 



V-9 

Thus the dual current density of Bosonic-field gravitomagnetic matter can be written as 

= y^"" - ^d-¥>*) = e'^'ja. (6.11) 

Note that here the index a refers to the one in the flat Minkowski spacetime. 

It is shown that in the above classical field equations of Fermionic and Bosonic field, gravitomagnetic monopole 
possesses the dual mass rather than the mass, namely, the concept of the mass is of no significance for the gravito- 
magnetic matter. Dual mass has its own gravitational nature and features different much from that of the mass. In 
this sense, matter may be classified into two categories, i.e., the regular matter and dual matter. The former has 
mass while the latter has dual mass. 

B. Dual current densities 

In accordance with the expression (6.7), the dual velocity (or current density) of dual mass can be defined as 
= e-i^Ua. To compare d^Uv — dvUn with duU^ — d^U^j, is physically interesting, which is valuable for us to define 
a new kind of the dual vierbein field. The calculation of d^Jj^ — d^U^ and dij,Ui, — dvU^ yields 

df,U^ - d^Uf, = [d^e^a - di,e^a - {e^^ad^ - e^adf^)] U". (6.12) 
Here we will adopt an alternative definition of the dual vierbein field e^a, i-s-., 

di^e^a - d^Bfj^a - {ei^adu - e^adij) = ^^ixv"^ [daefia - d/^eaa - {(^aadf) - C/Ja^a)] , (6.13) 

which is different from (6.1). Consequently, this definition will lead to the following duality relation 

di^eua - due^a - {e^i,adu - euadij) = -\^nv°''^ [dae^a - dfjCaa - {eaadff - e^ada)] ■ (6.14) 

Thus, it follows from Eq.(6.13) and (6.14) that Eq.(6.12) can be rewritten as dij,Uv — dvU^ = ^e^^"'^ [daUp — dpUa) 
and di^Ui, — di,U^ = — ^e^^"'' (daUp — dpUaj . So, the source tensors /C^,^ and H^v defined in Sec. V can also be 
expressed in terms of the dual vierbein field (or dual current density), i.e., there is a duality relation between tC^,^, 
Wjui/ and 'hLfxi,, 



The above discussion shows that if the vierbein field e^a is defined via Eq.(6.13), then there is no difference between 
the source tensor 5*^,^ and the dual source tensor 5*^,^. Here the dual source tensor S^u is the linear combination of 
IC^i, and Ti^v In a word, both the current density and the dual current density can be used to define the source 
tensor of gravitomagnetic matter in the gravitational field equation. 
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VII. THE FIVE-DIMENSIONAL CASE 



In order to find the potential connection between the gravitomagnetic charge and the magnetic charge, we will 
investigate the five-dimensional dynamics of gravitomagnetic matter. In the five-dimensional case, the dual Ricci 

tensor is 

-A _ 1^ CDE d fdgsc dgDB\ ,„ 

where A, B, C, D, E run over — 4, instead of — 3. In what follows Greek indices run over — 3. By using the 
cylinder condition {i.e., the derivative of the metric with respect to is vanishing), we can obtain 

^^4 = \e^ ^^^dr {dahx - dx9ai) + ^''^drdam- (7.2) 

By using the Kaluza mechanism, the five-dimensional metric qab is adopted as follows [20] 

where (j) and A^ denote a certain (constant) scalar field and the electromagnetic field, respectively. Thus Eqs.(7.2) 
can be rewritten as 

K'^ = + '"'drdaA,, = ^^^drT.x. (7.4) 

So according to the field equation (5.21), i.e., 7t^4 = -jp=T^4, the second equation in Eqs.(7.4) can be rewritten as 

= ^t^4, (7.5) 

which is actually the electromagnetic field equation of magnetic charge. This, therefore, means that the (/i4) com- 
ponent of the five-dimensional gravitational field equation of gravitomagnetic matter can be reduced to the electro- 
magnetic field equation of magnetic monopole via the Kaluza mechanism. Thus the five-dimensional dynamics of 
gravitomagnetic monopole provides us with a potentially unified way to treat the magnetic monopole and gravito- 
magnetic monopole. 



VIII. DUAL THEORY 



In electrodynamics, it is shown that the dynamics of magnetic charge is just the dual theory to that of electric charge 
(for a concise formalism, see Appendix 6 to this paper). This phenomenon arises also in gravity theory. Now we are 
therefore led to consider in detail the duality relationship between the dynamical theories of gravitomagnetic charge 
(dual mass) and gravitoelectric charge (mass). First we define a second-rank antisymmetric gravitational potential 
tensor hP'" as an alternative to the description of the gravitational field. Such antisymmetric gravitational potential 
tensor h^" may be so defined that the following expression 

e'*™'^ = cj{x) (/i'"^/i^^ K'^'W'' hyf'K"^) (8.1) 

is satisfied, where u{x) is a scalar function to be determined. 

Define a functional of h^" , say, Mr^v^i {h'^'^), which is assumed to contain the derivatives of h/^'^ up to second order 
{i.e., it contains deh^" , dgrjh'^'^)^. Based on this, we assume that the gravitational Lagrangian density is chosen to be 



It is assumed that the field equations associated with the gravitational field is often of second differential order. 
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By using the Euler-Lagrange equation, if the following equation is satisfied 



dC 



dC 



dL 



dh 



Act 



fin 



(8.2) 



(8.3) 



fiv 



then it is believed that the funetional A^^ 
keeping in mind Eq.(8.2), we can set 



Twi^/i (h^^) is related closely to the dual Riemann tensor TZt^vix, namely, 



TOJl/fJ. 



TUJl/fJ. ' 



(8.4) 



which results from the expression (4.2), |e'*™'"^xwi'M — 

The roles of Eq.(8.1), (8.3) and (8.4) is as follows: Eq.(8.1) can be used to determine the scalar function a{x); 
Eq.(8.3) is applied to the determination of the form of the functional Mtujvix {h'^'^); Eq.(8.4) can be employed to find 
the relationship between the antisymmetric gravitational potential tensor h^'^ and the symmetric metric tensor g^'' . 



Because of the identity, ^e'^^'^'^W 
this identity leads to i (G^'^ + G'^'^) 



1 ^Orbjy 



i (G^^ + G'^^), suggested in Sec. Ill, the substitution of Eq.(8.4) into 
'M^^,„^ {h^"^ . Thus Einstein's field equation is of the following form 



(rewritten in terms of the antisymmetric gravitational potential tensor h?^'^ rather than of the symmetric metric tensor 



with T^^ being the energy-momentum tensor of matter. In this sense, we have found the duality relationship be- 
tween the theories of mass and dual mass, and therefore developed a dual formalism (in terms of the antisymmetric 
gravitational potential tensor h^'^) which can express Einstein's equation within its framework. 



IX. CONCLUDING REMARKS 



Not all of the nonanalytic solutions of Einstein's vacuum equation belong to the contribution of gravitomagnetic 
charge. In an attempt to investigate the gravitational field caused by gravitomagnetic charge currents, we meet, 
however, with difficulties in selecting the "true" metric produced by gravitomagnetic charge from the solutions of the 
Bianchi identity {i.e., Einstein's vacuum equation). For this reason, one should have a dynamics of gravitomagnetic 
charge. To achieve this aim, this paper gives a concise presentation of the mathematical formalism of the field equation 
of gravitomagnetic matter. We derive the antisymmetric dual Einstein tensor from the variational principle with the 
dual curvature scalar, and construct an antisymmetric source tensor, which appears on the right-handed side of 
the gravitational field equation of gravitomagnetic matter. It may be believed that, as Lynden-Bell and Nouri-Zonoz 
stated [8] , the above-developcxl gravitational field equation might introduce a greater degree of physical understanding 
of NUT space (and Taub-NUT space). 

It is shown in this paper that the dual physical quantities such as the dual vierbein field (and hence dual metric), 
dual spin connection and dual currcmt density (and hence dual source tensor) constitute the dynamical framework 
of gravitomagnetic matter. It is also demonstrated that the topological dual charge needs its own gravitational 
field ecjuation. To clarify the relationship between the field equation of gravitomagnetic matter and Einstein's field 
equation is of physical essence. Einstein's field equation cannot govern independently the spacetime produced by 
gravitomagnetic matter. Instead, it can govern independently the gravitational field of regular matter only. The 
vaciuim equation of dual matter {i.e., without dual sources) is just the Bianchi identity of Einstein's equation, and 
similarly, Einstein's vacuum equation is just the Bianchi identity of the gravitational field equation of gravitomagnetic 
matter. The NUT solution is the one to both Einstein's equation and Eq.(5.18). 

It follows from the point of view of the classical field equation that matter may be classified into two categories, 
i.e., the gravitoelectric and gravitomagnetic matter. It is clear that the concept of mass is of no significance for the 
gravitomagnetic matter. In this sense, we can say that the gravitomagnetic matter possesses the topological dual 
mass [1,2]. 

Lynden-Bell and Nouri-Zonoz considered briefly the problem as to how general relativity must be modified if the 
gravitomagnetic monopole truly exists in the universe [8]. They thought that the space where the gravitomagnetic 
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monopolc is present may possess the unsymmetrical affinc connections [8]. Namely, their consideration means that 
the extra gravitational potentials should be introduced to allow for the gravitomagnetic charge. Indeed, the topolog- 
ical dual mass will introduce some additional new things in gravity theory, e.g., it will give rise to the topological 
singularities in spacctime, which is described by the nonanalytic metric functions. But there is no essential reason for 
cluing us on introducing the extra gravitational potentials. To deal with the dynamics of gravitomagnetic charge, the 
only thing we should do first is to introduce a new field equation {e.g., (5.17) and (5.18)), the vacuum case (without 
the sources) of which is just the Bianchi identity of Einstein's equation. 

Gravitomagnetic charge has some interesting relativistic quantum gravitational effects, e.g., the gravitational anti- 
Meissner effect, which may serve as an interpretation of the smallness of the observed cosmological constant. In 
accordance with quantum field theory, vacuum possesses infinite zero-point energy density due to the vacuum quantum 
fluctuations; whereas according to Einstein's theory of General Relativity, infinite vacuum energy density yields 
the divergent curvature of space-time, namely, the space-time of vacuum is extremely curved. Apparently it is in 
contradiction with the practical fact [21]. In the context of quantum field theory a cosmological constant corresponds 
to the energy density associated with the vacuum and then the divergent cosmological constant may result from the 
infinite energy density of vacuum quantum fluctuations. However, a diverse set of observations suggests that the 
universe possesses a nonzero but very small cosmological constant [21-24]. How can we give a natural interpretation 
for the above paradox? Here, provided that vacuum matter is perfect fluid, which leads to the formal similarities 
between the weak-gravity equation in perfect fluid and the London's electrodynamics of superconductivity, we suggest 
a potential explanation by using the cancelling mechanism via gravitational anti-Meissner effect: the gravitoelectric 
field (Newtonian field of gravity) produced by the gravitoelectric charge (mass) of the vacuum quantum fluctuations 
is exactly cancelled by the gravitoelectric field due to the induced current of the gravitomagnetic charge of the 
vacuum quantum fluctuations; the gravitomagnetic field produced by the gravitomagnetic charge (dual mass) of 
the vacuum quantum fluctuations is exactly cancelled by the gravitomagnetic field due to the induced current of 
the gravitoelectric charge (mass current) of the vacuum quantum fluctuations. Thus, at least in the framework of 
weak-field approximation, the extreme space-time curvature of vacuum caused by the large amount of the vacuum 
energy does not arise, and the gravitational effects of cosmological constant is eliminated by the contributions of the 
gravitomagnetic charge (dual mass). If gravitational anti-Meissner effect is of really physical significance, then it is 
necessary to apply this effect to the early universe where quantum and inflationary cosmologies dominate the evolution 
of the universe. Study of the geometric property in quantum regimes is an interesting and valuable direction. 

We think that there might exist the formation/creation mechanism of gravitomagnetic charge in the gravitational 
interaction, just as some prevalent theories provide the theoretical mechanism of the existence of magnetic monopole 
in various gauge interactions [25]. But so far it has to yield any deflnite theories of how the gravitomagnetic matter 
forms (we think that it may arise from the interactions associated with the Chern-Simons gauge field). The mag- 
netic monopole in the grand unified theory plays an important role in the early universe. The same may be said 
of gravitomagnetic monopole. If the gravitomagnetic monopoles truly exist in the universe, it will inevitably give 
rise to significant consequences and may also play an essential role in astrophysics and cosmology. In this paper we 
investigated the possibility of constructing a field equation that governs the gravitational field in the vicinity of this 
topological dual matter. We may conclude that both the dynamics of gravitomagnetic matter and other related topics 
associated with, for example, the creation mechanism of gravitomagnetic charge deserve further consideration. 
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APPENDICES 

Appendix 1. On the Levi-Civita tensor 

The Levi-Civita tensor e^^ap in the fiat Minkowski spacctime is so defined that it changes sign when any pair of 
indices are interchanged and £0123 = +1, e°^^^ = — 1. The Levi-Civita tensor eoi23 and e'^^^^ in the curved spacetime 
are defined to be eoi23 = -x/— 5^0123 and e°^^^ = --7==e°^^3, respectively, where g is the determinant of the matrix 5^^. 

One of the most important properties of the Levi-Civita tensor is that its covariant derivative is vanishing, which 
may be proved as follows: 

The covariant derivative of e^^gr is given 

^ A A A A 
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As an illustrative example, we will calculate eoi23:cr only. Here eoi23 = \/^.9^oi23 = \/— .9- It is easily obtained that 
—£0123 = ^\/—gg"^-^^gaf3- if t^^^^t = 0, l, 2, 3, respectively, then by using the antisymmetric property of e^iyOr, 



dx 

one can arrive at 



^crfAi/er 0<T^A123 — V-9^° Ocr> ^'^ jyCT^MAeT ^ icreOA23 — V^I"^ 1<ti 

9a^tlv\T —>■ 2cre01A3 = V~9^^ 2<t) ra'^fJ-^S^ ^ 3rTe012A = V^-^r^ 3<t • (^2) 

So, it follows from (Al) and (A2) that the covariant derivative of the Levi-Civita tensor €0123 is 

eoi23;. = I V^g'^" -^gais ' ^/^^^ = 0. (A3) 
An alternative derivation of the above result is given as follows: 

The transformation relation between the Levi-Civita tensors in the noninertial and local inertial frames of reference 
satisfies e^^0r = ^^%^^£abcd- So, the derivative of e^i^er reads 

def^^er _ d^"^ g^a g2^b g^c Q^d 

dx" dx^dx^^ dx" dx^ dx^ "''"^'^ dx^^ dx'^dx'^ dx^ dx'^ 

The insertion of = T\^^, ^ = r\,§^, ^ = and ^ = r\,§^ into (A4) yields 

de 9 

g^/ = CT/i^Ai/e-T + au^tJ^XOr + aO^iivXT + ^T-e^uex- (A5) 

It follows from Eq.(Al) and (A5) that enveT;<T = 0. Thus we have proven that the covariant derivative of the 
Levi-Civita tensor truly vanishes. 

In the four-dimensional spacetime, the following four identities (which are the inner products of two Levi-Civita 
tensors) are useful for investigating the dynamics of gravitomagnetic matter: 

^Xvad QlrtA 

^nva0^ — -o\g^, 
e^.a/je^""'' = -2! {g^gZ - g^gt) , 

e^uape^"^^ = - {g^Zgl + glgtgl + gW-gl - g^gigl - gW-ga - g^lal) ■ (A6) 



Appendix 2. Proving the first identity of Eq.(3.4) 
According to the fourth identity of Eq.(A6), one can arrive at 



^OtUJV A(T ^tOuJI/^ )LtA(T 



= - {g^'^g'^V^ + g'^'^g'^g'^^ + g'^^'g^^g'" - g^'^g^V - g'^'g^^g^'^ - g^'^g'V) ■ (A7) 

So, the further calculation yields 

= - {n<^>^ - \g>^<^n^ - {n''^ - \g^'^n^ . (as) 

If the Ricci tensor 7^^'^ and 7^'^'' are respectively defined to be 7e^^ = 7^''''^^, TZ''^'" = gxrH'^^''^, and the curvature 
scalar TZ is obtained via 1^^"^^ = —ga^T^"^ = —T^, then we can readily show that this identity 

1 ^eru..^^. = (G^'^ + G"^') (A9) 
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truly holds. 

Appendix 3. Showing that | ^/^e'^'^'^'' STZruvij.^'^ is a surface integral 
For the first we obtain the variation of TZruiv^ 

5 _ d 



+ ^jTy^^,. - ^JT^^,,. (Aio) 



For the second, we analyze the third term uv^^\,tij, on the right-handed side of Eq.(AlO). The following calculation 
is trivial 

= -e'^-r^jl (-^ + ^ + ^) = -e''-r\jT^,,,. (All) 

In the same manner, we can rewrite — e^™'T'*' ^^^^Fa.ti/ as follows 

_^^™.pA ^^^r^,.. = e''™^F^ ^^STrM- (A12) 
So, it follows from both Eq.(All) and (A12) that one can obtain 

^m™. ^^^p^^^ _ pA ^^^r^,,,) = e^™'^ (-r^ _<5r.,A^ + ^/r.,^,) , (A13) 

which will lead to some changes in the expression for t^'^'^" STZt^jv^- Thus the integrand e^'^'^'^dTZrui^ij. may be rewritten 
as 

/ g ^ ^ ^ 

(STr,^^)., - (^r,,„.).J . (A14) 
By using e^J'^'^ = e'^^'^'^ = 0, one can therefore arrive at 

ie^™'',57^™,^ = S'^^ = -^^(^S^) with = -e''™'^<5r,,,,. (A15) 

Hence, we have shown that | ^/^e^^'^'^'^STZruvudO, is truly a surface integral, i.e., 

I X = ^ (V^B^) dn, (AI6) 

which is a four-dimensional volume integral of the divergence (■^/^B''), and hence has no contribution to the 
derivation of the dual Einstein tensor. 

Appendix 4. The variation of the Levi-Civita tensor ^'^'^ 

With the help of the two expressions S^^^ = ^-^^ga^Sg"^ and dg^e = —gnageff^g"'^ , one can reach 
-^e/"" = S {g^ee'^''^) = Sg^ge'^^^ + g^eSe'^'^^ 



= { -g^.agefse'''^^ + ^g^egaoe'^^' ] Sg'^^. (A17) 
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Thus we obtain the variation of the Levi-Civita tensor e„ ^'^^ . 
Appendix 5. Mid-dual and two-side dual curvature tensors 

In Sec. Ill, we only considered the left-dual and right-dual Ricmann (and hence Ricci) ciirvaturc tensors, and 
showed that ^e^.'^^'^TZ^rXa = 0, ^Cj^ ^"^^T^-^xAo- = T^-va- In fact, there exist other kinds of dual Riemann tensors. See, 
for example, ^e^ '^^"'R-nrXa- It follows from the expression Te^^Aa = -^^i^,ra - -£^^n,r\ - r^r^./iA + r\^o,iJ,a- 
that 



T^HtXo- = -TlrtiXa + (^A^ct - d^dx) Q^r- (A18) 

Thus one can obtain 

^'^^n^rxa = {dxd, - d^dx) g^r- (A19) 

Note that here ie, {dxd„ - d„dx) g^^r = ^^"dxd^g^r = "^^^A (d^g^r - d^g^r) = -^H^r- So, a useful fact 
that no essential difference exists between ^e^ ^^'^'R-firXa and TZj^r is thus demonstrated. This, therefore, means that 
we need not take account of ^e^ ^^"Ti-nTXa in considering the dynamics of gravitomagnetic charge. 

In what follows, we will discuss the mid-dual and two-side dual curvature tensors. It will be shown that these two 
kinds of dual curvature tensors are either trivial or equivalent to the left-dual tensor. 

The mid-dual Riemann tensor under consideration is ^e^^^^'^TZ^xav, and the corresponding dual Ricci tensors can 
be written as 

5™ (^er..^"7^^A<x.) = 0, (^erc^'^/^Aa.) = 0, g^^ 0e.u.^"^/.Aa.) = nr., ••. (A20) 

It is readily verified that the mid-dual tensors in (A20) are either trivial or equivalent to the left-dual tensor. 

The two-side dual Riemann curvature tensor is defined to be ^e^^^'^TZxrua, the corresponding Ricci tensors of which 
are given as follows 

9""" (^S-^^'^T^a™-) = 0' (^e^/'^A^a) = -n.^., 5"" (^e^/^T^Ar...) = -Inur, ••. (A21) 

which are also trivial or equivalent to the left-dual tensor. 

For this reason, in this paper we will not consider further other kinds of dual tensors such as mid-dual and two-side 
dual curvature tensors. 



Appendix 6. The duality relation between the electromagnetic field equations of magnetic and electric charges 
The electromagnetic field equations of magnetic and electric charges are of the form (where and are the 
electric and magnetic charge current densities, respectively) 

d^J^>^''=r,, d^P^^=jl, (A22) 

with Tali = daAfi - dpAa, Tf"" = ie^''"^J^a/3, ^^"^ = -^e'"'°''^Taf3. The above electrodynamics is developed in the 
framework of the electromagnetic field In contrast, we can also propose a dual theory to the above one, where 
the formalism is constructed in the framework of the so-called dual electromagnetic field Af^, i.e., 

JT'*'' = a^i^ - a^i'*, Jr^i^^ = -^e^"'"'^Ta0. (A23) 

The above theories are dual to each other and, furthermore, equivalent to one another. This, therefore, implies that 
one can study the electrodynamics not only in the framework of A^ but also in the one of Jl^. In Sec. VIII, the 
similar dual phenomenon in gravity theory is discussed, where a second-rank antisymmetric tensor h^v is defined to 
describe the gravitational field. So Einstein's field equation can also be expressed inside the framework of h^v, just 
as the fact that the electrodynamics can be described within the framework of A^. 
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